arXiv: 1508.00186v2 [quant-ph] 30 May 2016 


Minimum Copies of Schrodinger’s Cat State in the 
Measurement of Multi-Photon Entanglement 

Yiping Lu 1 and Qing Zhao 1 * 

'School of Physics, Beijing Institute of Technology, Haidian District, Beijing 100081, PR. China 
"qzhaoyuping@bit.edu.cn 


ABSTRACT 


Multi-photon entanglement has been successfully studied by many theoretical and experimental groups. However, as the 
number of entangled photons increases, some problems are encountered, such as, the exponential increase of time necessary 
to prepare the same number of copies of entanglement states in experiment. In this paper, a new scheme is proposed based on 
Lagrange multiplier and Feedback, which cuts down the required number of copies of Schrodinger’s Cat state in multi-photon 
experiment and keeps the standard deviation error of fidelity unchanged. It reduces five percent of the measuring time of 
eight-photon Schrodinger’s Cat state compared with the existing, and also guarantees the low error in fidelity. Additionally, the 
same approach has also been applied to the simulation of ten-photon entanglement. It reduces the twenty two percent of the 
copies of ten-photon Schrodinger’s Cat state compared with the uniform distribution, and the distribution of optimized copies 
gives better fidelity estimation than the uniform distribution of copies in the case of the present the ten-photon Schrodinger’s 
Cat state. 


Introduction 

From fundamental tests of quantum mechanics' to quantum teleportation, quantum key distribution, and quantum communi¬ 
cations,®] quantum entanglement has wide applications in different areas. Recently, a single photon has been recognized to 
teleport multiple degrees of freedom simultaneously, which includes spin and orbital angular momentum. 5 The Greenberger- 
Horne-Zeilinger (GHZ) states created in experiment have been obtained by combining the momentum and polarization. 6 — 
Several experiments have been performed to validate multi-photon entanglement. 6 8 10 12 An indispensable tool is the en¬ 
tanglement witness for certification of entanglement in these related experiments. 13 15 Generally, the expectation value of 
entanglement witness can be evaluated by fidelity.® Precise estimation of fidelity requires many identical copies of the prepared 
state. 17 On the other hand, the coincidence count rate of multi-photon entangled states decreases exponentially with a linear 
increase in the number of entanglement photons, which is generated by the phenomenon of a nonlinear process of parametric 
down-conversion in BBO. 18 22 Hence, collecting sufficient copies of multi-photon entanglement state costs much longer time, 
such as, in eight-photon entanglement it takes 170 hours to produce the sufficient copies of eight-photon Schrodinger’s cat (SC) 
state (See the label of Fig.3 of RefPj). Up to now, a study on ten-photon entanglement or more is inaccessible in experiment 
since the coincidence count rate of ten-photon entanglement state is lower than 9 counts per hour. 7 It needs nearly three months 
to prepare sufficient copies, for example, 110, of the ten-photon SC state to certify entanglement according to the current 
technology. (See Appendix “Preparation of Ten-photon SC state”). 

Discrimination of a quantum state by adaptive process is developed recently. Adaptive process is to split the conventional 
measurement into several pieces and to choose the current measurement based on the results of previous measurements. 
The standard of selection is chosen to be minimizing the probability of error. 23 The probability is estimated by the known 
information. Generally the adaptive process is spl it into two steps. The first step is to get rough information. The second step is 
to rectify it and get a precise density matrix. 24 25 

In this paper, an efficient method is developed to reduce the number of copies of an unknown state to certificate entanglement 
in the multi-photon experiment. Specifically, the conventional measurement is split into several steps. For each step, the least 
number and the optimal distribution of identical copies of the unknown state on different measurement settings are calculated 
by the proposed model. The measurement result from previous steps provides the value of parameters for future steps. In 
our model, the unknown state is supposed to be a pure SC state or SC state in the presence of noise. Since the entanglement 
validation of SC state is through fidelityjS the optimization introduces fidelity as a criterion. When fidelity is greater than 0.5, 
the experimentally prepared state is certified to be entangled. 13 The target of optimization is to search for the minimum number 
of copies of the unknown state that can confirm the error gap of fidelity belonging to a small area; therefore the fidelity interval 
can be estimated and the minimum number of copies of state is obtained. 





Minimum copies of multi-photon Schrodinger’s Cat state 

The experimental n-qubit SC state is denoted by a 2" x 2" density matrix p exp . Its fidelity with the pure state |SC) is defined as 
F exp (\SC)) = ( SC\p exp \SC) = Tr( Pexp \SC){SC\). (1) 

To calculate the F exp (| SC )), liq.ljTJ) can be written as 
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where (w) is the expectation of entanglement witnessE^^D Hence, F exp (\SC)) can be calculated by evaluating (w). In Eq.([3jl, 
|SC)(SC| is decomposed into the form 
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where M kn / n = cos(kn/n)o x + sin(^^/w)(J J . lcl HI See Appendix “Entanglement Witness” for more details. 

The n-qubit SC state requires at least n + 1 settings to calculate fidelity (see Observation 1 Based on Eq.(j5]), the 

standard deviation of fidelity is deduced. 
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In Eq.([6j», tj is the total number of copies of n-qubit entanglement state that projected into the jth measurement setting, 
j = 1,2, • • • , n + 1. Its value equals the sum of accumulated n-fold coincidence counts in all different bases of the /th setting. 
Here accidental coincidence count is ignored since it is almost zero when n is large. The P\ is equal to the summation of two 
relative frequencies. One is the case that all qubits are projected into horizontal polarizations (|//)(//|)®", the other is the case 
that all qubits are projected into vertical polarizations (|V)(V|)®". Here, the meaning of relative frequency is the ratio between 
the number of copies of a state projected into a base and the number of copies of the state measured in all the bases belonging 
to this setting. Similarly, the Pj is the linear combination of relative frequencies of different basis in the y'th setting. It should be 
kept in mind that measurement setting means a group of complete basis that copies of a state are projected in the same period 
of time and relative frequencies being gained simultaneously. The details are shown in the appendix “Standard Deviation of 
Fidelity”. 

We intend to apply fewer copies of the unknown state to estimate fidelity with same accuracy. Let £q denote the given upper 
bound of standard deviation of fidelity since the number of copies of a state required relies on it. Our objective is to use as few 
copies of the state as possible and, at the same time, to narrow down the fidelity to a small interval. Therefore, the following 
model is proposed: for «-qubit SC state, we have 
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It should be noted that tj obtained by solving Eq.(|7]) is sufficiently large, since the larger tj is, the higher the probability for 
the result of Eq. 0 to hold, which will be discussed in section of “Characteristics of optimization of the successful probabilities”. 
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Based on numerical results, the solution of Eq.(|7jl has large tj in most cases and the probability for the above model is nearly 1. 
Following is the analytical solution of Eq.([7j) obtained. Let Sq = £, then 
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The process to get the Eq.([8]) can be found in “Theoretical derivation of minimum copies of multi-photon Schrodinger’s Cat 
state” 


Results 

Direct estimation of fidelity for experimental eight-photon SC state and simulated ten-photon SC state 

The advantage of our method over the existing approaches can be demonstrated by the experiment of eight-photon entanglement. 
When sufficient copies of eight-photon SC state in the presence of noise (Psphotons) are projected into different settings 
in experiment, fidelity is calculated by total accumulated coincidence counts on different basis and then the eight-photon 
entanglement can be verified. 7 In this section, our method is to change the number of copies of prepared SC state ( Psphotons ) 
measured in different settings. The results show that total copies of prepared SC state can be saved; while, fidelity precision 
remains same. 

To apply Eq.(j8]) for certificating of eight-photon entanglement in our model, some parameters need to be given, such as 
n = 8. Let the prepared eight-photon SC state in experiment be psphotons , which is a SC state mixed with noise. In order to 
have optimized results compared with the experiment. The error bound of fidelity, £o, is set at 0.016, which is exactly the 
same value as the one used in experiment. This number can be found in the second to last paragraph of Ref. 7 According 
to entanglement witness, an eight-photon SC state requires nine settings to determine fidelity uniquely, see equation (2) in 
RefPI Let | H) represent horizontal polarization and |V) represent vertical polarization. And define |+)=(|//) + e l8 \V)) /\/2, 
|—)=(|//) — e ,e \V))/y/2. Then the nine measurement settings are defined as 


V« = {(l+,0)(+,0|f 8 ,(l+,0)(+,eir(|-,0)(-,e|),---, 

(|-,0)(-,0|)® 8 ,0=O}, 

= {(|+,©)(+,0|)® 8 ,(H-,e)<-F,0|)® 7 (|-,0)<-,0|),---, 

(|-,0)<- 0|f 8 ,0 = 7T/8}, 

SsphotonsA = {(| + ,0)( + ,0|)® 8 ,(| + ,0)(+,0|)® 7 (|-,0)(-,0|),--- , 

(|-,0)<- 0|)® 8 ,0=27t/8}, 

58 Pfc5 = {(|+,0)(+,0|)® 8 ,(|+,0)(+,0|)® 7 (|-,0)<-,0|),---, 

(I-,0)(- 0|)® 8 ,0 = 37t/8}, 

58pfc6 = {(|+,0)(+,0|)® 8 ,(|+,0)(+,0|)® 7 (|-,0)<-,0|),---, 

(I-, 0)(-,0|)® 8 ,0=47T/8}, 

58 P / !O r o „ s ,7 = {(|+,0)(+ ) 0|r 8 ,(l+,0)(+,0|)® 7 (|-,0)(-,0|),---, 

(|-,0)(-,0|)® 8 ,0 = 5tt/8}, 

V« = {(l+,0)(+,e|f 8 ,(l+,0)(+,0|) 07 (|-,e)(-,e|),---, 

(|-, 0 )(-, 0|)® 8 ,0 = 6 tt / 8 }, 

V« = {(l+,0)(+,eir 8 ,(l+,0)(+,0|) 07 (|-,0)(-,e|),---, 

(|-,0)(-,0|)® 8 ,0 =7tt/8}. 

(9) 

Furthermore, 1305 copies of psphotons are prepared in total in experiment. 7 Notice that 1305 is not directly given in the 
Ref., 7 but is used to draw the graphes, calculate fidelity in Ref. 7 and provided by the author of that paper. The number can also 
be roughly calculated by the copies of p&photons P er h° ur and the total hours spent. That is9x40 + 9x25 + 9x 15x7 = 1544, 
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Setting 

Experiment 

Uniformity 

Optimization 

S 8 photons,! 

352 

145 

415 

S 8 photons,2 

200 

145 

106 

S 8 photons,3 

107 

145 

103 

S 8 photons,4 

100 

145 

106 

S 8 photons,5 

110 

145 

103 

S 8 photons, 6 

111 

145 

108 

S 8 photons,! 

106 

145 

101 

$8 photons, 8 

116 

145 

108 

S 8 photons,9 

103 

145 

103 

Summation 

1305 

1305 

1253 


Table 1 . Distribution of copies of psphotons under different case. The first row represents the number of copies of Psphotons that 
is the summation of the number of accumulated coincidence counts projected into all the bases of the first measurement setting. 
The following eight rows represent the number of copies of p% p hotons that projected into the other eight settings corresponding 
to 0 = 0, ?r/8, • • ■ ,77 t/ 8. The last row is the total cost of copies of Psphotons i n the experiment of Ref., 7 uniform distribution 
and our optimization. The second column is the cost of number of copies of p^photons in the experiment for different setting. 
The third column represents uniform distribution of copies of Psphotons in eac h setting. The last column represents the copies of 
PSphotons obtained in the optimization. 


in which the coincidence counting rate can be found in the 11th paragraph of Ref. 7 and the hours spent for different settings 
can be found in the label of Figure 3 of Ref. 7 Accidental coincidence count is very small in eight-photon experiment, therefore 
it is neglected. A numerical test is performed. Firstly, a set of copies of a quantum state measured at various settings is defined 
as “distribution of copies”. Three different distributions (experimentally applied distribution of copies of psphotons, optimal 
distribution of copies of the state, which is obtained from Eq. and uniformity distribution of copies of the state) are 
considered separately, and compared with each other. The number of copies of Psphotons for each case is listed in Table[l] The 
first column is the tag of setting. The optimal distribution of the copies of Psphotons calculated by Eq.([8j» is listed in the last 
column. Obviously, the total number of copies of psphotons required is cut down to 1253, thus 52 copies of Psphotons (about 5 
percent) are saved compared with experiment. Since coincidence count rate is only nearly nine (8.88) per hour (in the 11th 
paragraph of Ref. 7 ), around 5.9 hours could be reduced in the experiment while the same precision of fidelity can be hold. 

For each case, fidelity can be calculated from new relative frequencies obtained by simulating the experimental process in 
computer according to the real precise relative frequencies in different settings. In simulation, the real relative frequency is 
calculated according to Born’s rule. It also requires density matrix to be known in this rule. Fortunately, the density matrix 
of Psphotons can t> e obtained by experimental data and phaselift approach, which will be given in detail in “optimization of 
multi-qubit experimental and simulated data via density matrices”. Since summation of the real relative frequency of different 
bases in a same setting is equal to one, the interval between 0 and 1 is divided into 2 8 sub-intervals and the range for each of the 
sub-interval is equal to the value of the corresponding relative frequency. And a random number between 0 and 1 is produced 
with the equal probability for each value between 0 and 1. And the interval it lies in is found and the number of event for 
this interval is added to one. After producing random numbers with the number of copies of p% p hotons for the setting, different 
interval gets a different number of event. Then relative frequencies can be calculated. After that simulated fidelity is obtained. 
We also divide the fidelity range from 0 to 1 into 50 equal portions. Event number is added to one when the calculated fidelity 
belongs to the corresponding interval. All three situations (experiment, optimization and uniformity) are all repeated for 550 
times separately, which means 550 fidelities are calculated. The number of events per interval is accumulated and observed, as 
shown in FigurejT] Figure [l^ shows that: when all 1305 copies of psphotons are applied, the experimental results give better 
estimation of fidelity than the uniform distribution since the height of the outline for uniform distribution on vertical axis 
direction is lower than the experimental one. The outlines for experiment and optimization are also described, which almost 
coincide with each other. However, optimization only costs 1253 copies of the Psphotom^ which is smaller than the 1305 copies 
of P&photons required by the experiment, as shown in Figure [lj). The Figure [lj: demonstrates that optimization is also better than 
the uniform distribution. 

At present there is no way to create enough copies of ten-photon SC state to certify entanglement in experiment. Numerical 
test is produced to estimate fidelity based on a computer created density matrix Piophotom that its fidelity with pure ten-photon 
SC state is 0.8414. It is carried out in the situation of uniform distribution in each setting (100 copies of ten-photon SC state for 
each setting) and optimization. The process is same as eight-photon entanglement. Both cases are all repeated for 100 times 
separately, then the distributions of fidelities are obtained, as shown in Figure [2] It is observed that 22.45 % copies of simulated 
ten-photon SC state Piopiwtons can be saved according to Eq.([8j> compared with uniform distribution on each setting. 
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Obviously, the optimization yields a better estimation of fidelity with limited copies of state available. The scheme given 
here is useful in certifying the multi-qubit entanglement state and can be generalized to any state by changing the form of 

constraint of Eq. (13- 


Optimization of multi-qubit experimental and simulated data via density matrices 

In addition to the direct estimation of fidelity, we also estimate a density matrix first by phaselift, 27 and then calculate fidelity. 

The model for calculation of the density matrix is constructed based upon the procedures given in Refs.® - 55 the noise case 
is applied, which is in 36 and, 27 

Minimize £ |Tr(pM M . v ) -f py \ 

M-v 

subjectto p>0,Tr(p) = l, (10) 


where p is density matrix, M p y is positive operator valued measure (POVM) in the p-th bases of the v-th setting, f py is the 
relative frequency in the p-th bases of the v-th setting. 

The quantum state tomography for three, four and eight-photon entanglement is conducted. When corresponding experi¬ 
mental frequencies f pv are put in Eq.( 10 1 , the density matrix is calculated out. Our objective is to use the least copies of an 
unknown state to obtain a density matrix close to the real one. The real density matrix p exp is approximately obtained with the 
use of a large number of copies of the state prepared in experiment. Then, p exp is applied to gain new frequencies according to 
Born’s rule through the simulation of experiment process on computer. These frequencies are applied to obtain the density 
matrix p re . Then many density matrices p re are obtained under different number of copies of the state and compared with the 
p exp achieved in the experiment. 

Several examples are given below in the following. The density matrix of the three-photon SC state ( Pexp ^qubits ) is obtained 
by using the experimental data and construction method summarized by Eq.([T0|) with Pauli measurement, as shown in Figure 
[3ji. The two large elements on the diagonal of the density matrix p exp -} qu hii\ are equal to 0.50188 for \HHH)(HHH\ and 
0.38419 for |VVV) (VVV |. And the real parts of two main elements on the anti-diagonal are both 0.37238 on | HHH) (VVV| and 
\VVV)(HHH\. The imaginary parts are quite small, so are not drawn. The density matrix of four-photon SC state ( Pexp-4qubits ) 
is also obtained by using experimental data and phaselift, as shown in Figure [ 3 J 5 . The result is obtained by Pauli measurement 
and Eq.([To|). Two large elements on the diagonal of the density matrix equal to 0.50637 for | HHHH) (HHHH\ and 0.36161 
for \VVVV) (VVVV|. The real parts of elements on the anti-diagonal are 0.35944 on | HHHH) iyVVV | and \VVVV)(HHHH\. 
Since the p eX p- 3 q ubits (Figure [3^) and p exp _4 qu bits (Figure[3j)) have very small noise and the purity is high, density matrix of 
three-qubit ( p3 qu bits ) (Figure |3p. Figure |3jl) with much more noise is created for the following simulation. Two large elements 
on the diagonal of the density matrix are equal to 0.3716 for \HHH)(HHH\ and 0.3412 for |yyy)(yyy|. And the real parts 
of two main elements on the anti-diagonal are both 0.3504 on \HHH) (V VV | and \VVV)(HHH\. The corresponding imaginary 
part is nearly approach to zero, so is not shown. The density matrix of the eight-photon system ( p& p i w tom ) is a ls° drawn in 
Figure [4] Obviously, only the real part of elements in four corners of the density matrix are larger than 0.2; other elements 
are much less than it, which is the characteristic of SC state. Besides, the imaginary part is too small; therefore, it is not 
drawn. By the density matrix of Figure [3ji, we reconstructed the p exp _ 3 qubits under different number of pauli measurement. The 
reconstructed density matrix is p e st. Figure [5]exhibits fidelities and Mean Square Error (MSE) when a different number of 
POVM is applied. When sampling number of POVM achieves around 45, fidelity is in a stable value (around 0.7) and the 
corresponding MSE is near 0. 

Similarly, the following four measurement settings are defined for three qubits measurement. 


Sw*s,i ={m(H\f\m(H\f 2 (\V)(V\),---,(\V}(V\) & } > 

S3 q ubi t s, 2 = {(\+,e)(+,e\r\(\+,e)(+,e\r 2 (\-,e)(-,e\),---, 

(\-,o)(-,o\f 3 ,o = o}, 

53, H M il 3 = {(l+,0)(+,0|)® 3 ,(l+,e)(+,0ir 2 (|-,0)(-,0|),---, 

(|-,e)(-,e|f 3 , e = 7T/3}, 

5 3 ^,4 = {(|+,0>(+,0|)® 3 ,(|+,0)(+,0|)® 2 (|-,0)(-,0|)v, 

(|-,e)(-,e|f 3 ,e = 27r/3}. 

(ii) 

Two measurement settings shown in, 37 which is to project the state into S^ qu bit s ,i and S-i qu bits,h can also be applied to certify 
entanglement. For the three-qubit SC state, the density matrix ( P 3 qu bits ) are measured in four settings or two settings, respectively, 
and then the fidelity is estimated, as shown in Figure [6] In the figure, it shows the distribution of number of event of fidelity 
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between a target pure SC state and estimated one when 10000 copies of three-photon SC state piqubits are measured. The fidelity 
between Piqubits and pure three-qubit SC state is 0.7068. We use Random to represent the distribution of 500_ 1000_5000-3500, 
which represents the number of copies of piqubits prepared in four settings, such as 500 is prepared in the Si qu bit s , \, 1000 is 
for the setting of Sj qu bi ts , 2 , 5000 for Siqubits ,and 3500 for S^qubitsA- Optimization means 3630-2570_2670_ 1130, in which 
3630 copies of Piqubits is projected into the setting of Sj qu biis. 1 , 2570 copies of p^qubiis is projected into Si qu bit s , 2 , 2670 is for the 
setting of S^qubitsp^ 1130 is for the setting of S^qubitsA- “Uniformity ” means 2500_2500_ 2500-2500, which means all four 
settings are projected with the same number of copies of Piqubits (2500). “Two setting” means 5000-5000, which represents 
the number of copy of p-iqubits prepared in two settings, such as 5000 is for S^ qu bits, l, an d 5000 for S^qubUsp- Both two-setting 
distribution and the optimized distribution of copies of Piqubits (3630-2570-2670-1130) gives the best estimation of fidelity, 
while the randomized distribution (500-1000-5000-3500) gives the worst estimation. Christ mentioned that a bias exists for 
fidelity estimation when the semi-definite constraint is added to the maximum likelihood approach, and this bias is based on 
density matrix. 38 Here phaselift is applied, and there is no obvious bias for fidelity estimation when the number of copies of 
Piqubits approaches 10000, as shown in Figure [6] However, there is an obvious bias when the number of copies of p-iqubits drops 
to 1000 and the number for the setting of Siq U bits .2 is switched into the setting of Si qu bitsA in two settings case, as exhibited in 
Figure [7] 

Fidelity estimation is also compared and analyzed in different initial conditions, such as the number of copies of Piq U bit s , as 
shown in Figure[8] which presents that 10000 copies of piqubits provide much better estimation of fidelity than 200 copies of the 
state. Furthermore, optimization always gives better estimation of fidelity than uniform distribution of copies of piqubits- 


Optimization of the number of copies via experimental feedback 

Let us note that Pj in Eq.(|8j should be known before calculating tj, and there is no way to obtain the precise value of Pj without 
knowing density matrix via Bom’s mle or without experimental measurement. In “Direct estimation of fidelity for experimental 
eight-photon SC state and simulated ten-photon SC state”, we prior estimate a density matrix based on the preparation of 
copies experimentally. In “Optimization of multi-qubit experimental and simulated data via density matrices”, we estimate a 
precise density matrix by phaselift. Here, we show to calculate them through the experiment itself. Take eight-photon SC state 
experiment as an example. In experiment, pure eight-photon SC state is the target state that needs to be prepared. It can be 
taken as a priori to approximately decide Pj, such as P\ is a value near to Tr((|//)(Z/|)® 8 |SC)(SC|) + Tr((|V)(V |)® 8 |SC)(SCj), 
Pjis near to a value given by Eq.< 301, y" = 2,3,--- ,n+ 1. However, the experimentally prepared state is not pure \SC) and it 


takes too long time to precisely estimate Pj, an optimization procedure is proposed based upon the experiment. It divides the 
process of measurement into a few steps. Instead of measuring one setting for a prolonged time to obtain the frequencies within 
small error margins, and then continue to measure the next setting for the same time, and so on. We divided this total long time 
into several intervals, and changed the order of measurements. The order is to measure all the required settings one by one for a 
much shorter time, then based on the measurement results; the Pj can be estimated roughly. After that, the extra number of 
copies of a quantum state that needs to be prepared and measured for each setting can be given by Eq.([8]l by inputting the rough 
Pj. Afterwards, more copies of the quantum state are prepared and measured according to the tj given. Later, more precise 
frequencies can be obtained and this process can be repeated until the final precision for fidelity is reached. Figure |9]shows the 
measurement order when the process is conducted only twice. We simulated this process in computer, it only costs a very short 
time, as shown in Figure [TO] 

To be specific, the main process is as following. Firstly, we introduce a superscript to represent the number of steps in 
optimization. The superscript / of a parameter represents the parameter applied in the Z-th step, i.e. e 1 represents the value £ 
used in Eq.([8]» for the first round of measurement. At the beginning of fidelity estimation, £ 1 is set to a large number, such as 
0.01, and all of Pj are originally set to Pj = Pj = 1/2, j = 1,2,3, • • • ,n + 1, (Pj can also be chosen according to target pure 
state |SC), such as P) can be a value near Tr((|//)(Z/|) 0 "|5C)(5C|) +Tr((|y)(V|)®"|5C) (SC\) so that the suitable solution tj 
can be obtained by solving Eq.(|8j). The experiment is performed according to the fj copies of the quantum state. When all tj 
copies of the prepared quantum state is projected into a measurement setting, Pj can be obtained. After that, input Pj instead of 
Pj , and have £ 1 become smaller; consequently, the tj can be obtained. Then copies of the state with the number of tj — tj are 
projected into the jth measurement setting in the second round of experiment, so on and so forth. Measurement is ended when 
E ller is sufficiently small, obvious, 


l 1 3 

£ > £ > fc > 


> £ 


> £ 


( 12 ) 


Extra time is needed for optimization; however it is much shorter compared with the time required for the preparation of 
copies of multi-photon entanglement state, as shown in Figure [10] The iteration makes the experiment to have more pauses 
between different settings during the measurement process. Mostly, switching settings cost more time, generally is about 3 or 4 
times of the switching time in conventional measurement. Anyhow the time required for optimization is much shorter than that 
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spent on the preparation of the copies of multi-photon entanglement state. Generally, switches and optimizations only cost less 
than two minutes, while the coincidence count rate of eight-photon entanglement state is too low that it costs several hours 
to produce just enough copies of the state for only one setting. Therefore, the time in calculation and switching time can be 
neglected compared to the preparation of copies of multi-photon entanglement state. 

In the following, a specific example is given. Numerical simulation applies a four-qubit SC state mixed with gaussian noise. 
The density matrix is PAqubits , which fulfills trace equal to one and semi-definite condition. The fidelity between PAqubits and 
pure SC state is 0.9374. In simulation, parameters are chosen as e 1 = 0.01, £ 2 = 0.001, £ 3 = 0.0001, £ 4 = 0.00001. The five 
measurement settings are required and listed as follows: 

SAqubitsp ={(\H)(H\r\(\H){H\r\\V)(V\),--- ,(|F)<V|)® 4 }, 

VM i ,2 = {(l+,0)(+,0|f 4 ,(l+,e)(+,0|f 3 (|-,0)(-,0|),---, 

(|-,e)(-,e|f 4 , e = 7T/4}, 

VM i ,3 = {(l+,0)(+,0|f 4 ,(l+,e)(+,0|f 3 (|-,0)(-,0|),---, 

(|-,0)(-,0|f 4 ,0=27r/4}, 

(|-,0)(-,0|f 4 ,0 = 37T/4}, 

v^,5={(i+,0)<+,0ir,(i+,0)<+,0ir 3 (h,0}<-,0i)r--, 

(|-,0)(-,0|f 4 ,0 = ^} 

(13) 

All the initial numbers of copies of PAqubits f° r eac h setting are set at 5. Other initial parameters for 5 measurement settings 
are set at: Pi = 1 /2, Pi = 1 /2, P3 = 1 /2, P4 = 1 /2, P5 = 1 /2 , respectively. In Figure[TT[ the PAqubits is taken as the test matrix. 
Its horizontal axis represents the number of iteration, which means the number of round of measurement. The corresponding 
point is the average number of extra copies of the state PAqubits that needs to be projected into each setting for the next round of 
measurement. The curve connects the number of required copies of PAqubits f° r eac h same setting. The error bar is one standard 
deviation, which is obtained by repeating the optimization program for 100 times. When the iteration ends, the Pj is listed in 
Table [2] 


Pi 

0.9339 

p 2 

0.0316 

p 3 

0.9491 

Pa 

0.0325 

P 5 

0.9474 


Table 2. The final Pj when iteration ends. 


We define the ratio between the e 1 applied in the current round of measurement and the £ / 1 applied in the previous round 
of measurement as the ratio of £, and let the ratio of £ for different round of measurement be equal to each other, that is 
£ 2 /E 1 = £ 3 /£~ = ■■■ = £ / / / £ / 1 = ■ ■ ■. The different values are tested to search for the best ratio that costs the least number of 
copies of a state. Figure [12] shows the number of copies of state randomly created at different ratios of £ when it equals to a 
value between 0.05 and 0.9. It is observed that it requires most number of copies when the ratio is 1/2, it increases at a wave 
type before this value, and decreases after this value. Specifically, the following procedure is conducted. Initially, the number 
of copies of a randomly created state is an integer between 4 to 7 for each setting. Pi, P 2 , • • •, P 5 are all given a value randomly 
created between 0.25 and 0.75 and £ is 0.01. Then, a fixed ratio of £, such as 0.1, is applied. It means £ is set to 0.001 in the 
second round of measurement, 0.0001 in the third round of measurement and so on. Eq.([ 8 ]i is applied to calculate the needed 
copies (tj) of the state for each round of measurement. Then tj is summed up to gain the total number of the copies of the state 
in current round of measurement The iteration stops when £ is smaller than 0.0003. The minimum number of copies of 

the state can be found by repeating the above steps by changing the ratio of £. It is found that the value near 0.15 for the ratio 
requires the least number of copies of a randomly created state, which is 178. The total iteration number for each setting for 
most created states is about 3 or 4 to achieve the 0.0003 of finial £. It is also noticed that the required number of copies is even 
less when the ratio of £ approaches 0.9. However, it is not suitable to apply the large value for the ratio of £ since too less 
copies of a state may lead to our model hold with a low probability, which will be discussed in the following section. 
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Discussion 


Characteristics of optimization of the successful probabilities 

It is noticed that the summation of number of copies of an unknown state projected into the same setting must be larger than a 


certain value, since a large value can confirm the model to hold with probability nearly one. 

In the above sections, the minimum number of copies of an unknown state is obtained and its fidelity belongs to the 
interval with the certain high probability. Hoeffding’s inequality is a mathematical way to describe the probability. It states 
that the sample average X = J^X/t of t independent, not essentially identical distributed, bounded random variables with 
Prob[Xj £ [ai.bi}\ = 1 for i = 1,2, • • • ,t satisfies 


Prob[X-E[X] < -h] 

< exp[-2rh 2 /£(&,■ - a,) 2 ] 


(14) 


for all h > 0, where X, is a variable, a,- is the lower bound, /;, is the upper bound, t is the number of samples, E[X] denoting the 
mean value of X, h is the definite value that equals to the maximum deviation from expectation/ 9 40 

Now this inequality is applied to multi-photon entanglement certificate experiment. The measurement response of a single 
copy of an unknown state is taken as the value of a single random variable. Since photon detectors can only give the feedback, 
0 or 1, this leads to /?, = 1 and a, = 0. X corresponds to a relative frequency, which is denoted to be fj, j is to distinguish 
different measurement settings. Hence the expectation E[X\ corresponds to the probability p r The total copy of a state for the 
j,/, setting is represented by tj instead of t. By replacing all of them, we obtain 


Prob\fj-pj < —hj\ < exp[-2tjhj\, 


(15) 


where hj is the deviation from true probability pj. It means 


Prob[fj £ (pj-hj.pj + hj)\ X \-2exp[-2t j h 2 j]. 

Then 


(16) 


Prob [1 ~fj£ (1 — Pj “ hj. 1 - P j+ hj)\ 
^ 1 — 2exp[—2tjh 2 }. 


(17) 


Therefore, 


Prob \fj — (1 -fj) £ (2(pj-hj)-l,2(pj+hj)-l)\ 


^ 1 —2exp[—2tjh 2 ]. 
Let fj (1 — fj) be pj. Hence, 


(18) 


Prob [Pj £ (2 ( Pj - hj) - 1,2 ( Pj + hj) - 1)] 
/ 1 — 2 exp[—2tjh 2 \ ,j = 2,3, • • • ,n +1. 


(19) 


Let 2(pj — hj) — 1 be P- and 2(pj + hj) — 1 be Pj~, and based on Eq. (34 1 .kj£ [,fe- ,k+], where 



( 20 ) 



( 21 ) 



( 22 ) 



(23) 
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then from Eq.([8]), one has 


tj = 


( 24 ) 



(25) 

Therefore the f, e [t■ ".// ] in Eq.([8j» when the holding probability of model Eq.<j7j» is considered. 

Obviously, hj has the impact on f ( + and tf . The larger hj is, the larger the gap between tj and /■ . Large hj and t, from 


Eq.( 15 i are needed to keep results with high probability. However, large t, costs too much experimental time. Large hj 


may introduce too large a gap between r and T , which may then lead to the wrong number of copies of an unknown state. 
Therefore, it requires to choose suitable hj and tj. 

By comparing the optimization results with the experiment, it is found that only 986 copies of pg qu bit s are used compared 
with the 1305 copies of pg qu bits in eight photon experiment, which specifies that 24 percent copies of eight-photon SC state 
(pgqubits) can be saved. Specifically, hj is chosen to be 0.2 for all j. According to joint probability, JE pj is calculated, in which 
Pj is the successful probability for each setting. For eight photon measurement, j = 1,2, • • • ,9, the final probability is 0.9972 
for experiment after 1305 copies of pg qu bits are measured. We observed same probability is obtained when 110 copies of pg ql ,bits 
for each setting are used and all /z ; s’ are chosen as 0.2. 

In the above analysis, we assume Hoeffding’s inequality describes the probability precisely. In the following, numerical 
simulation is produced to confirm the above mathematical tool is true. The density matrix ( pg qu bits ) is calculated from 
experimental frequencies, and new relative frequencies are obtained under a certain number of copies of pgqubits in a random 
simulation of experimental process that gets the relative frequency by computer. 1 — Pi is the summation of relative frequency 
that all the qubits projected into horizontal polarization and the relative frequency that all the qubits projected into vertical 
polarization. The real value of 1 — Pi is 0.8068 when the number of copy is sufficiently large. When failing probability is set 
less than 0.0001, Figure [j~3|shows how 1 — Pi behaves under different number of copy of pg q „bits- In the figure, red circle and 
blue triangle are drawn according to Hoeffding’s inequality, 1 — Pi can be estimated much more precisely with an increasement 
of prepared copies of the state pg ql ,bit s - It is observed that all the numerical simulated points lie in the region between upper 
bound and lower bound. Therefore, Hoeffding’s inequality can be applied to describe the Pj in multi-photon entanglement. 


Extension of the optimization of the number of copies of a state to quantum-state tomog¬ 
raphy 

The surprising thing brought to us by the optimization in the fidelity estimation, is to extend it to quantum state tomography. 
The optimization model for tomography is constructed as follows: Let po be a d x d density matrix of real experimental created 
and p be the estimated density matrix via limited copies of po. For n qubit state tomography, we build 

Minimize ^ T v 

V—1 

subject to | |p — po| \ F < £o, (26) 

where T v represents the number of copies of po of the vth setting, v can distinguish different measurement settings, n s is total 
number of measurement setting. The solution of Eq.(p6| is 

See “Theoretical derivation of minimum number of copies of a state in quantum-state tomography” for the details to get the 
solution of Eq.(|26|) 


Conclusions 

We proposed an optimal approach that assists to find the minimum distribution of copies of a state that is sufficient to certify the 
entanglement of the state by fidelity. The main purpose is to facilitate an experiment to obtain better measurement strategy for 
fidelity estimations, for example, by changing the ratio of the number of copies of the state in different settings. To estimate 
fidelity directly from fewer copies of SC state (1253 copies), with optimized distribution, almost the same distribution of fidelity 
as the experimental one (1305) can be obtained. It not only saves time, about five percent of measurement time (6 hours) is 
saved, but also small error of fidelity. Additionally, the distribution on the number of copies of ten-photon SC state is also 
simulated, and 22.45% of copies of the ten-photon SC state are saved, which further highlights the superiority of this scheme, 
and reveals that the optimized distribution of copies of a state in different settings gives better estimation of the fidelity than 
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uniform distribution of copies of a state in all settings. Fidelity can also be estimated by the reconstructed density matrix. It is 
observed that the optimized distribution provides the best estimation of the true state, the uniform distribution provides a worse 
estimation, while randomized distribution provides the worst estimation. With the increase of the number of copies of the state 
the differences between different distributions (uniform distribution and optimized distribution) become much smaller. Besides 
the state with high similarity with SC state, this approach can also be extended into other states in parallel. And the scheme is 
extendable to tomography when the MSE between the estimated density matrix and real density matrix is limited to a fixed 
value. 


Preparation of ten-photon SC state 

From second paragraph of Reference,® the count rate of eight-photon event is about 2.8 x 10 5 Hz. Accidental coincidence 
counts can be neglected for eight-fold entanglement. Therefore two-photon event count rate is \Z2.8 x 10 S Hz . Detecting 
ten-photon entanglement requires totally five independent pairs of entangled photons to present at the same time, so the ten- 
photon coincidence event scales as (v/2.8 x 10 ~ 5 x 3600) 5 = 0.0568 per hour. For ten-photon entanglement, 11 measurement 
settings are required according to the entanglement witness of SC state. If only 10 copies of ten-photon SC state are 
prepared and measured in one setting, then 110 copies of SC state are required. Therefore, the corresponding time is 
(110/0.0568 )hours = 1.9366 x 10 3 hours = %0.69Y1 days ~ 3 months. 


Entanglement witness 

To calculate F exp , each term in the decomposition of \SC) (SCj has to be measured to determine its expectation value. For an 
eight-qubit SC state, n = 8, the expectation values of all the terms on the right hand of Eq. <0 should be calculated. Specifically 
the total accumulated coincident counts on the i-th base is defined as n,s such as n\ copies of ps p hotom- with all qubits are 
projected into horizontal polarization | H). «256 copies of p% p hotons with all qubits are projected into vertical polarization \V). 
Relative frequencies on (|//)(//|)® 8 or (|F)(V|)® 8 can be calculated by «i/(L?=i n i) °r «256/(L?=i «,)■ To get the expectation 
value of the third term of Eq.Q, we have 

= [—Tr PexpM^J s + TrpexpMf^ - Tr pexpMfn/S 

+Tr PexpMf£/ s — TTpexpMf*j s + Tr PexpM®*i£ 

-Trp exp M^+Tr Pexp Mf] 

= iE (-!)*<</ 8 >> (27) 

k= 1 

in which (M^ 8 g ) represents the expectation of the operator Mj^ 8 g . The estimation of expectation value of the operator M^. 8 g = 
(|+, 0) (+, 0\ — | — ,9)(— ,0|)® 8 is equivalent all the expectations of various combinations of |+,0)(+,0| and | —, 0)(— ,0|, 
due to 

= (l+,0)(+,eh|-,e)(-,0|)® 8 

= (|+,0)(+,e|f 8 -(|+,0)(+,e|)® 7 (|-,0)(-,0|) 

+.... + (|-,0)(-,0|)® 8 . (28) 


There are 256 terms in all for a fixed 0. When the number of copies of state that projected into different combinations of bases 
|+, 0) and |—, 0) are collected, the copy numbers corresponding to Af^ 8 , g are calculated from Eq.(28 1 . Thus, the (M^ 8 , g ) can 

be evaluated . 10 13 From these measurements, the expectations of different terms appearing in the decomposition of the SC state 
entanglement witness are obtained. 


Standard deviation of fidelity 

The error is calculated from Poisson distribution in Fig.2 and Fig.3 in Ref . 10 for each term of Eq.(|5]>. 

Based on the experimental data of Ref ., 7 the all eightfold coincidences are mainly projected into (| H) (H |)® 8 or (|V)(V |)® 8 in 
S9,photons, I setting. When the state is projected into the setting of horizontal or vertical polarization, the copy of Ps c/ ubits projected 
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into (|//)(/7|)® 8 is 148, the number of copies of ps qu bits is projected into (|V)(V |)® 8 is 136. The summation of total number is 
68 for the case that eight qubits are projected into other bases in Ss p hotons,i ■ Therefore, the ratio Pi between the number that 
projected into (\H)(H\)®* or (|y)(y|)® 8 and the total number is (148+136)/(136+148+68)=(148+136)/352=284/352=0.8068, 


the ratio (1 — Pi) between the copy that some qubits are projected into the horizontal polarization | H) and some are projected 
into vertical polarization |V) and the total copy of p^photons f° r this setting is 1-0.8068=0.1932. The smaller value in P\ and 
1 — Pi is defined as P\. Similarly, according to Eq.(281. A smaller value between (M®. 8 -^g) and 1 — is chosen 


as Pj, j = 2,3, - - - ,9. Therefore, P 2 = 40/200 = 0.2, P 3 = 20/107 = 0.1869, P 4 = 20/100 = 0.2, P 5 = 21/110 = 0.1909, 
P 6 = 23/111 = 0.2072, P 7 = 19/106 = 0.1792, P 8 = 24/116 = 0.2069, P 9 = 20/103 = 0.1942, in which the largest ratio is 
0.2072. 


When the number of copy of state in the whole measurement time are large and the relative frequency that the copy of state 
projected into a base or several bases in a setting is close to 0, Poisson distribution can be approximated by binomial distribution 
(Page 291 of Ref. 41 ). Notice that the Poisson distribution here is not for the entangled photons created in BBO in time scale but 
the distribution of number of copies of state on different measurement basis satisfied. The binomial distribution is a special case 
of the Poisson binomial distribution, which is a sum of tj independent non-identical Bernoulli trials. 42 In our optimization model, 
binomial distribution is applied since the number of copy of state is quite large. Let P\ represent the probability that the copy of 
state is projected into \H)(H\® n or |V)(V|®" in the setting of Si, where Si = {\H)(H\® n , \H)(H\® n ~ l |V)(V |, • • • , |V)(V|®"}. 
And let Pi denote the ratio that the state collapses to other bases in Si, hence Pi = 1 — Pi. According to Fig.2 and Fig.3 of 
Ref., 10 a value in Pi or P\ is close to 1 and the other is close to 0 when t\ is much larger than 20. It satisfies the condition 
that Poisson binomial distribution can be approximately replaced by binomial distribution. Since the variance of the binomial 
distribution is fiPi(l — Pi) (Page 277 of Ref.®), then variance of number of events that SC state collapses to \H)®” or \V)® n is 
also the same value since fiPi(l —Pi) = /1Pj (1 —Pi). Therefore the standard deviation is ^/t\P\{\ —Pi). Besides, the Pi is 
defined as the ratio between the number of copies of state detected on a (\H)(H\)® n or (|V) (V|)® n basis and the total copies of 
state in Si. Therefore the standard deviation for the relative frequency is y/t\P\[\ —P\)/t \, which is equal to \JP\ (1 - P\ )/t\. 

From Eq.ijTJ and Eq.([5]», 


Fexp{\SC)) 

= Tr( Pexp \SC)(SC\) 


— Tr ^ Pcxp 

Tr ^ Pexp 
= Tr < p exp 


* + Tr < Pexp 


\ ( (\H)(H\f" + (|V)(y|r + - £(—1 )*M®'/, 

Z \ n k= 1 

i((|//)(//|)®" + (|y)(y|) 

(|//)(//|r + (|y)(y|) 

\®n 


\®n\ 


\®n 


— y (-i 

7 „ i—i y ' kit n 

k=l 


k=l 


L rr S Pexp =-(-1) 


= Lr{p, w [(|H)(H|)* + (|l')(V|)®"]} + £ P(-l)*Tr{p„, [m“ J } 

k= 1 ^ 

= ^Tr{p„ /> [(|ff><ff|)*" + (|y)(y|)®*]} 

£ ^-(-i) A ‘Tr{p avp [(|+,0)(+,e|)®"-(|+,0)(+,e|)®"“ 1 (|-,0)(-,e|) + --- + (|-,0)<-,0|)®"]} 

:=1 LU 

^Tr{p eX p [(|7/>(//|)®" + (|y)(y|)® n ]} 

£^Tr{p CTp [(|+,0)(+,0|)®" + (|+,0)(+,0|)®' I - 2 ®(|-,0)(-,e|)® 2 + . ■• + (!-, 0)<- 01)®”]} 

k= 1 m 

£^Tr{p e , p [(|+,0)(+,e|)®' l - 1 (|-,0>(-,e|) + (|+,0){+,0|)®' I - 2 (|-,0)<-,0|)®(|+,0)(+,0 

k =1 


■ 2 
n 

+ E 

k= 1 

= -Tr 
2 


l)+- + (|-ex- 
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Define 


A = Tr { p exp [ ( \H) (H | f” + (| V) (V | H } , 

pj =Tr{p« p [(i+,0)(+,0ir+(i+,0)(+,0ir“ 2 ®(i-,0)(-,0|f 2 +. ••+(!-, 0)(-,eir]}, 

i-A = Tr{p e , p [(| + ,0)(+,0|r“ 1 (|-,0)(- I e|) + (|+,e)( + ,0|r“ 2 (|-,0)(-,0|)®(|+,0)(+,0|) 

+...+(i-,0}(-,0ir- i ®(i+ ! 0)(+ 1 0i)]}}. oo) 

Then Eq.(p9|) can be rewritten as 


1 

2 


n 


* + I 

k= 1 


(- 1 )* 
2 n 


A+l — 


n 


E 


(-i)* 

2 n 


(i-flt+i) 


1 

2 

1 

2 


^E^.-E^ 


*=1 


*=1 

«+i r_iv'-i «+i 

^E^-E 

7=2 7=2 


2 n 


(31) 


Here Eq.(|28|> is applied when n = 8 and k + 1 is denoted as j in last second step. 
According to the previous analysis, the standard deviation of Pj is 


/A(i-A) 


, 7 = 1,2,---,9. 


(32) 


Further considering the formula of combined standard uncertainty , 43 the standard deviation of fidelity can be derived. We 
use A F exp to represent it. Therefore, 


AA 


exp 


\ 


dF exp / A(l-A) 
dPi V ?i 


n +1 

E 

7=2 


dF exp Pj(l-Pj) 


dPj 


\ 


1A(1-A) , 1 y A(! - A) 


4 ?! 


7=2 


(33) 


Theoretical derivation of Minimum copies of multi-photon Schrodinger’s Cat state 

Let 

A = ^A(l - A), 

k J = pl 2 P j( l - P j)’j = 2 ’ 3 ’--- >«+ L 

e = el (34) 

then the optimization problem is equivalent to 


n +1 

Minimize tl jt2) ... ltn+1 £ ?; 

7=1 

«+! £ . 

subjectto ^ — < e,? 7 - > 0,kj > 0,j = 1,2,-• • ,n+1 (35) 

7=1 

where e and kj are positive real constants, n is a positive integer, tj is a variable and also a positive integer. In order to solve 
Eq.([35]> easily, all the variables, including the number of copies of state, tj, are considered as a real. The optimized number of 
copies of state is then rounded off to the smallest integer greater than the final real tj. 
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_i-1 k ■ 

The solution of the optimization problem is assumed to satisfy L , ,-=i jr = £■ If the optimal solution is not on the boundary, it 

J 1 J 

means yf < £. Appropriate reduction in the number of tj can be made, while the inequality (X/i! < e) is still satisfied. 


This is contradictory with the target function “Minimize tj”, therefore the optimal solution must exist on the bound. 


7=1 J 


The detailed process of how to find the analytical solution of Eq.(35 i will be shown below. 
Now let 


tj — 


1 


(36) 


By substituting Eq.(|36| into Eq.(|35]l, we obtain 

n+1 j 

Minimize ^ 


j—i x j 


subject to 
Y J kj x j = e- 


(37) 


The Lagrange multiplier method is applied to solve the problem. Since the target is the minimization of £"jjj jr, Eq.(|37 


leads to 




(38) 


To find its minimum, partial derivative for each xj is expressed as 

dL 1 , „ „ 

jX ~ 


(39) 


From Eq.(39 1 , the following equation is obtained, 

^ 

The constraint of Eq.([37]> is 
Y J k j x j = £- 


From Eq.(40 1 and Eq.(41 1 , A is given by 

2 


A = — 




(40) 


(41) 


(42) 


By substituting Eq.(|42]> into Eq.(40 1 , we arrive at 


Xi = 


J y/kj(SjtWk])' 

Rewriting Eq.(43 i using Eq.(f36|), we have 


(43) 




y/kj&r = Wkj) 


(44) 


By substituting Eq.(|34[i into Eq.(|44|. Eq.([ 8 ]i is obtained. The optimal results of Eq.([ 8 ]) can be compared with experiment 
when the same coefficient P\ and Pj s are substituted. 
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Theoretical derivation of minimum number of copies of a state in quantum-state tomog¬ 
raphy 

According to Bom’s rule, 

P[i,v — Tr(Mjj :V po)) 
fn,v=Tr(M^ v p), 

H = 1,2,--■ ,d,v = 1,2,--- ,n s , (45) 


where p distinguishes different measurement operators in the same setting, d is the dimension of density matrix and P^y is the 
probability when po is measured by operator M LL V . Namely, p € C d/d , then 

Po = 

ft.v 

P = ^//i.v44u.y (46) 

ft.V 

Therefore one has 


p-pO = J^(fn,v-Pn,v) M n,v 
M,v 

Consider /W u v is orthogonal to each other, then 

IIp-Po||f 

= [Tr((p-p 0 )(p-po )*)] 1/2 


(47) 


^C/ju,v ~ P/r,v) 2 Tr(M |UiV M^ v ) 

M.v 


1/2 




T H,v 


Tr (M^yM^y) 


1/2 


(48) 


In the last step of Eq.(|48|>, standard deviation of binomial distribution is applied, see the first paragraph of “Standard 
deviation of fidelity” in appendix for detail. When the measurement operator M^y belongs to the same setting v, they have 
the identical number of copies T^ v of po, i.e.7/ v = Tiy = ■ ■■ = Tj v . We denote them to be T v . The model of Eq. (26 1 is 
equivalent to 


n s 

Minimize ^ T v 

V— 1 


subject to 



/ai,v(1 — /m.v) 


1 1/2 


Tr {M^yM^y) 


< So- 


(49) 


It is easy to find that the target of Eq.(|49|> is similar to Eq.([7]» except the larger required number of settings and different 
coefficients. 

Let L£ = i[/ M)V (1 ~ fn,v) Tr(M^ jV M^ v )] be k v and £ = £q, then the model has the following form 


Minimize ^ '//■ 

V—1 

" ' A: 

subject to ^ — < £,T V > 0,£ v > 0,v = 1,2, ••• ,n s . (50) 

v^\ Lv 
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( 51 ) 


Obviously, it has a similar form to Eq.(|35j); therefore the solution is the same as that of Eq.([ 8 ]> 

%/MLv=i Vkj) 


T v = 


If M^y is non orthogonal with each other, then 

IIP —Polk 

= {Tr[(p —p 0 )(p Po )*]} 1/2 


E (fn,v-Pp,v)(f ll 'y -^y) Tr (M^ v M*, y ) 


/ / 
,V 


/ I 

,v 


\UA'-Uy) 


l py 


\ 


f , ,(l -f, ,) 
J /i ,v r ,v > 

r / / 

U .v 


Tr(M M!V M% v ,) 


1/2 


Applying the similar substitutions as used in the orthogonal case, we have 


Its 

Minimize E T v 
v—1 

«j k 

subject to E 


v,v 


vy=t^VV 

r v >0,70 >0,k v i >o, 


< e, 


V ’ v,v 
V 1 5 2, • • • ,72^,V 1 j2, • • • 

Substitute q v = 1/7V, the constraint in the optimization becomes 

A. / 


v y =l 

n s 

= E 


v,v =1 
n s 


< E *Vv'(3V + 4 v ')/ 2 


V.V =1 


<E 

P=1 

< e. 


E v + E 


.V =1 


V —1 


q P A 


(52) 


(53) 


(54) 


Therefore, the non orthogonal case has a similar result with the orthogonal one Eq.(50 1 except the coefficient is different. 
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Figure 1 . a: The number of events versus fidelities for both experimental distribution and uniform distribution, b: The 
outlines for experiment and optimization, which almost coincide with each other. Optimization only costs 1253 copies of 
Psphotons • which is a smaller number than the 1305 copies of psphotons required by the experiment, c: The number of events 
versus fidelities for both optimization distribution and uniform distribution. 



Fidelity 


Figure 2. The distribution of events’ number of fidelity of ten-photon entanglement state ( P\o p hotons )■ The blue vertical line 
represents the fidelity between the simulated state P\o p hotons and pure ten-photon SC state, which equals to 0.8414. Different 
lines are used to connect adjacent points. Black line represents the optimized fidelity distribution and red line represents 
uniform distribution. It is observed that optimized distribution has more events accumulated near the real fidelity at 0.8414 than 
uniform distribution. 
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Figure 3. Density matrices iPexp-iqubits, Pexp-Aqubits > P:\qubii.P- Different colors are applied to represent the value of elements 
of density matrices, a: The real part of experimental density matrix of three-qubit ( Pexp-iqubits )■ b: The real part of 
experimental density matrix of four-qubit ( p eX p-4qubits )■ c: Real part of three-photon density matrix ( p3 qu bits) by random 
creation, d: Imaginary part of three-photon density matrix ( p3 qu bits) by random creation. 
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Figure 4. The real part of experimental density matrix of eight-photon SC state ( Psphotons )■ Different colors are applied to 
represent the values of the elements of density matrix. 
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Figure 5. Fidelities and MSEs’ of state p esl under different number of POVMs. p est is the estimated density matrix from 
Pexpqubits- The fidelity is calculated by Fidelity es t,pure = Tr(p est \SC) (SC\). The MSE (Mean square error) is calculated by 

PISE T f(p es i Pexp-3qubits)^ {Pest Pexp-3qubits)\- 
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Figure 6. The distribution of fidelities between the target pure SC state and the estimated states under different number of 
copy distribution on different settings. “ Random ” means the number of event of fidelity when distribution of copy of p^qubits 
goes “500_1000_5000_3500”, in which the 500 copies of state p^qubits are projected into Sj qu bit s . i; 1000 copies of pj,q U bits is 
projected into S^qubitsX 5000 copies of p^qubits is projected into the basis of set of S^qubiisp and 3500 copies of p^qubits is 
projected into S^qubitsA- Similarly, “Optimization” represents the event number of fidelity when distribution of copies of p^qubits 
is “3630-2570_2670_ 1130”, in which 3630 copies of state p^qubits are projected into S^qubitsX 2570 copies of p^qubits is 
projected into S 3qubits x, 2670 copies of p 3? „ to , s is projected into S 3qub it s , 3 a nd 1130 copies of p^qubits is projected into S 3qubUsA . 
“Uniformity” means the distribution is “2500-2500-2500-2500”, which represents all is equal to 2500 for the number of 
copies of state p^quba.s that projected into S 3qub i lSi u S 3qub i ts , 2 , Siqubitsp, and S 3qub i, S A ■ “Two Setting” represents the both equals to 
5000 for the copies of p3 qu bits that projected into S^qubitsA and ^'Squbits.i ■ The range of fidelity is split into 250 intervals on 
average between 0 and 1 to compare the event number. Different number of events that fidelity lies into a certain interval is 
gained, such as, if the calculated fidelity is 0.005, then it belongs to the interval between 0.004 and 0.008, the number of events 
belong to this interval is added to 1, and so forth. Fidelity is estimated for 550 times in all four situations. Black squares 
represent the number of events accumulated in each interval for the case of “500-1000-5000-3500”; Red circles represent for 
the case of “3630-2570-2670-1130”; Blue triangles represent for the case of “2500-2500-2500-2500” and pink triangles 
represent for the case of “5000-5000”. The fidelity between pure three-qubit SC state and p^qubits is 0.7068. The points are 
connected by the lines with the same color of the points. It is observed that red circles and pink triangles have the most events 
near this value. Therefore both the “Optimization” (“3630-2570-2670-1130”) and “Two setting” (“5000-5000”), perform 
better for the estimation of fidelity. While the black squares for the random one gives the worst estimation. 
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Figure 7. The distributions of fidelities between target pure SC state or experiment p eX p-3qubits and estimations obtained from 
the 1000 copies of Pi qub i ts - The cyan vertical line at 0.7068 represents the fidelity between p 2qub us and pure three-qubit SC 
state. Green vertical line at 0.82274 represents fidelity between experiment p eX p- 3 qubits and P ure SC state. In the figure, “a” 
represents the fidelity between the estimation and experiment p exp _ 3 q U bits ■ “b” represents the fidelity between the estimation 
and target pure state. “499_1_1_499” represents the distribution of fidelity when 499 copies of p^qubits are projected into the 
bases of S 3qublts _, ; 1 copy of p^ qublts is projected into the bases of S 3qub u S: 2 , 1 copy of P 3 qllb its is projected into S iqubits ^ and 499 
copies of p^qubits are projected into S^qubits ,4. The number of events of fidelities between estimations and target pure state in this 
case is denoted by pink triangle, pink dashed line is applied to connect them. The number of events of fidelities between the 
estimation and experiment in this case is denoted by blue triangle and blue dashed line is applied to connect them. 

“363 .257 267 113” represents fidelity distribution obtained from optimization distribution on four settings applied by 1000 
copies of experimental three-qubit SC state (p eX p_ 3 q ubits) built by Pauli measurements. “363 257 267 113” represents 363 
copies of p 2qubits are projected into S 3qubitS} 1 , 257 copies of pj, qubits is projected into S 3qubitSt 2 , 267 copies of pi qubi t s is projected 
into S 2qub i ts3 and 113 copies of p^qubits is projected into S^qubits, 4 - Red circles are applied to denote the events of fidelities 
between the estimations and experiment. Red line is applied to connect them. Black squares represent the number of events of 
fidelity between the estimation and target pure state. Black line is applied to connect them. It is observed that 
“363~257_267_ 113_a” performs better for fidelity estimation than “499-1 l_499_cz”. “363_257_267_ 113_£>” performs 
better than “499_ 1_ 1_499_7>”. Therefore, optimization distribution of copies of p 2qub i ts performs better than the performance 
of two settings. 
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Figure 8. The distributions of number for events of fidelities between estimations and target pure SC state when the fidelities 
are obtained from the density matrices that are constructed by 200, 1000 and 10000 copies of Piqubits- The black down triangle 
connected by black solid line represents the number of events of fidelities when they are calculated by optimization distribution 
of total 200 copies of p^qubits on all the four settings. Purple star connected by dashed line represents the number of events of 
fidelities when it is calculated by uniform distribution of 200 copies of P 3 qu bits ° n all the four settings. Blue up triangle 
connected by solid line represents the number of events of fidelities when it is calculated by optimization distribution of 1000 
copies of piqubits on all the four settings. Magenta diamond connected by dashed line represents the number of events of 
fidelities when it is calculated by uniform distribution of 1000 copies of p^ qu bits on all the four settings. Olive square connected 
by solid line represents the number of events of fidelities when it is calculated by optimization distribution of 10000 copies of 
Plqubits ° n a ll the four settings. Red circle connected by short dashed line represents the number of events of fidelities when it is 
calculated by uniform distribution of 10000 copies of p^q U bit s on all the four settings. Optimization distribution of 200 copies of 
Plqubits is “73 51 53 23”, which represents the distribution of fidelity when 73 copies of state P;\ q „i n , s are projected into 
S-iqubits. I; 51 copies of pi qubits is projected into S 3qub u s , 2 , 53 copies of is projected into S iquhlls 3 and 23 copies of p^qubas 

is projected into the bases of S^qubitsp- Similarly, uniform distribution (“50 50 50 50”), (“250 250 250 250”), 
(“2500_2500_2500-2500”), optimization distribution ( “363_257_267_ 113”) and (“3630-2570-2670-1130”) all follow 
the same rule as “73-51-53-23”. Namely, the first number is the number of copy of p 3qu bits that projected into the setting of 
S'Squbiis. 1 ; the second number is the number of copy of pj qu bits that projected into the S 3qu bit s ,z\ the third number is the number 
of copy of pq, qu bits that projected into S 3qu bit s p and the last number is the copies of p 3qu bus that is projected into the bases of 
S^qubitsp- 10000 copies of p 3qu bits gi ve much smaller error or standard deviation of fidelity than 200. Optimization always gives 
more centralized estimation of fidelity than uniform distribution under the same number of copies of p 3qu bits- 
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Figure 9. Comparison between traditional and optimization measurement of three-qubit state. The color of line segment 
represents the different setting. The length of line segment represents the time for the corresponding measurement. Traditional 
measurement order is to finish the measurement of each setting one by one, as shown by the Pastel yellow area. The 
optimization measurement is iterated twice as shown by the light blue area. 



Figure 10. The required time of optimization for different number of copies of a state. The red circle represents running time 
of eight-photon optimization. Black square represents running time of four-qubit optimization. It is observed that the total time 
to calculate Eq.([ 8 } and simulate the experiment is less than 100 seconds for both cases. Therefore compared with several hours 
spend to prepare copies of eight-photon state in experiment. It can be negelected. 
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Figure 11 . The change of optimization number of copy of p^qubits of fidelity estimation corresponding to different iteration 
numbers with different settings: k/ 4, 2n/4, 3zr/4, and K. “Setting for H/V” represents the number of copy of P 4 q U bi ts projected 
into Scubas, i • “Setting for 6 = 7T/4” represents for S 4 qu bu s p, “Setting for 6 = 2n/4” represents for S 4 qu bi ts 3 , “Setting for 
0 = 3k/4 ” represents for S^qubitsA, “Setting for 0 = n” represents for S 4 qu bu s , 5 - Error bar represents one standard deviation. 



Figure 12. The required copies of a state for different ratio of e. For each ratio of e, numerical test is conducted for 10 times. 
Black square is used to represent the number of copies of the state. Error bar is mean standard deviation. It is observed that the 
number of copies of state rises at wave type when the ratio of e increases but no larger than 1/2 and decreases with the ratio of 
£ when it is larger than 1 / 2 . 
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Figure 13 . The change of 1 — Pi in the setting of .S’i corresponding to different number of copies of eight-photon SC 
stat e(p squ bits)- Black square represents the numerical simulation, red circle represents the theoretical lower bound and upper 
triangle represents the theoretical upper bound. Numerical simulation is repeated for 10 times for each number of copies. It is 
observed that all the simulated points lie in the region that consists of point that is larger than the lower bound and smaller than 
the upper bound. 
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